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ABSTRACT 



Natural convection liquid cooling of simulated electronic compo- 
nents in a vertical channel was investigated. The test surface 
contained a single column of eight rectangular, protruding heated ele- 
ments, each simulating a 20-pin dual-in-line package. Temperature 
measurements and flow visualization were performed for a number of 
power dissipation levels and channel widths. Collectively, this infor- 
mation was used in interpreting the flow and transport characteristics. 
A correlation to predict the heat transfer rates was developed based 
on the component surface temperatures. Optimum channel widths 
were determined from these surface temperature measurements for 
the range of power levels investigated. Temperature distributions in 
the fluid were measured using a traversing thermocouple probe. 
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I. INTRODUCTION 



A. STATEMENT OF PROBLEM 

The continuing trend toward microminiaturization of electronic 
components has made thermal design critical to achieving higher 
packaging densities. The ever-increasing heat fluxes must be removed 
while keeping the device junction temperatures typically below 100°C 
[Ref. 1]. The failure of microelectronic chips increases exponentially 
as the junction temperatures increase [Ref. 2]. 

Microchips act as discrete sources of thermal energy. Effective 
removal of the generated energy is of major import in determining 
how they might be used without overheating the elements themselves 
or the surrounding package. Thermal control is one of the critical 
considerations in the design of any electronic cooling system. The 
thermal designer’s task is further complicated by the presence of 
arrays of electronic components placed in proximity to each other, 
requiring detailed considerations of transport interactions between 
neighboring elements. 

The design of packaged electronic components must consider all 
anticipated modes of heat transfer. A complete thermal analysis will 
account for all three mechanisms of thermal energy transport: 

• conduction 

• convection 

• radiation 
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The areas of concern include transferring heat, via conduction, from 
internal heat sources to the external surface of the package and the 
subsequent removal of heat from the external surface by a cooling 
medium. Also, the coolant temperature must be maintained at pre- 
scribed levels for a given environment. 

The present work focuses on the heat removal from electronic 
packages through convective cooling. This convection may, in general, 
be forced or natural; the medium liquid or gas. Because of the higher 
density and thermal conductivity of liquids compared to air, heat 
transfer coefficients an order of magnitude or more larger than for air 
are attainable. Due to the significantly high heat transfer rates achiev- 
able, cooling of electronic components by direct liquid immersion has 
been identified as an extremely effective thermal control technique 
[Ref. 2], 

Convection heat transfer rate from a surface at temperature T s of 
area A to surrounding fluid environment Tjnf is given by Newton’s law 
of cooling: 



q = Ah(T s - T^) 

This is a definition of h, the convective heat transfer coefficient, more 
than a phenomenological law. Evaluation of h depends upon heat 
source geometry, orientation, fluid properties, and fluid velocity. 

Natural convection results when fluid density gradients exist in 
the presence of a body force field such as gravity. Density gradients in 
the fluid commonly occur due to temperature gradients. Advantages of 
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employing natural convection in electronic cooling include high relia- 
bility, low cost, reduced noise, and a simpler design. Despite these 
potential benefits, only a limited number of studies are available on 
liquid immersion natural convection cooling of electronic components. 

Most existing natural convection studies, both experimental and 
analytical, have considered heating surfaces that were semi-infinite. 
Electronic component cooling configurations typically employ arrays of 
small, discrete, flush, or protruding heat sources. The available heat 
transfer correlations are inadequate in predicting transport in these 
geometries. Convective heat transfer from an array of microchips is 
strongly dependent on package shape and array density. Indeed, in 
real systems, individual elements or modules often vary significantly in 
size, shape, and in the manner in which they are arranged on the 
substrate. 

The geometric complexity is further compounded by the nonuni- 
form nature of heat dissipation from the array. Individual component 
dissipation is a function of its utilization. Variation in power dissipa- 
tion from component to component on the same board and from board 
to board within the system usually exists. The orientation of the 
boards on which the elements are mounted is yet another variable. 

Often, these semi-regular arrays face the back of an adjacent 
board. The resulting vertical channel has one surface which is rela- 
tively smooth and the other one covered by an arrangement of large, 
heat-dissipating elements. Within the channel, local convective heat 
transfer is driven by two phenomena, a local buoyancy force due to the 
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heated components and a forced convection effect arising from the 
flow history, over the entire channel length, including both heated and 
unheated regions. A simultaneous consideration of these mechanisms 
is essential in obtaining predictive capability for transport rates in 
these systems. 

The complexity of the problem described demands an approach 
which builds on the fundamental physical processes. Efficient cooling 
cannot be attained without understanding the heat transfer from each 
package by a determination of its flow and thermal fields. The overall 
goal is to develop a methodology based on sound physical principles 
which will allow the prediction of individual component temperatures, 
by superposition, in a nonuniformly heated array of components. 

B. IMMERSION COOLING— ANALYTICAL AND EXPERIMENTAL 

STUDIES 

Reference 3 provides an excellent literature review of previous 
work, both analytical and experimental, in the field of liquid immer- 
sion cooling of electronic equipment. Additionally, very relevant work 
has been performed by Ortega and Moffat [Refs. 4-6]. They have 
experimentally examined natural convection air cooling of simulated 
electronic components protruding from a vertical surface, with and 
without a shrouding wall. Heat transfer rates significantly higher than 
for a smooth vertical plate were found. They interpreted this as an 
indication of turbulent free convection. When a shroud was in place, 
they found plate-averaged heat transfer coefficients 40 to 50 percent 
higher than for a smooth, parallel plane channel. 
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Jaluria [Ref. 7] performed two-dimensional computations for 
natural convection air flow over a vertical surface with multiple wide, 
flush-mounted heaters. The study treated the flow as a boundary layer 
problem. Rajakumar and Johnson [Ref. 8] reported computations of 
free convection heat transfer using rectangular strip heating surfaces. 
Their finite element approach showed significant stagnant areas 
between the blocks. 

C. OBJECTIVE 

An experimental study was undertaken to to determine the 
hydrodynamic and thermal characteristics of buoyancy-induced flow 
over several heated protruding components mounted on a vertical 
surface. 

A single column array was chosen because it provided the sim- 
plest geometry featuring discretized heat dissipation. The effects of a 
neighboring card were simulated by a movable vertical shroud plate. 
Specific objectives of this study were: 

• to visualize the transient and steady-state, natural convection flow 
within an interrupted channel, for a range of component power 
dissipation levels and channel widths: 

• to measure component temperatures for various power inputs and 
shroud spacings to develop an appropriate nondimensional cor- 
relation to predict the heat transfer rates for this geometry: and 

• to measure temperature distributions in the adjacent fluid, using a 
traversing thermocouple probe to assist in evaluating the hydro- 
dynamic and thermal nature of the flow. 
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This study was a follow-on to work described in Reference 3 and 
is part of a continuing effort to study natural convection liquid immer- 
sion cooling of electronic equipment. 
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II. EX PER I MENT 



A. GENERAL DESIGN CONSIDERATIONS 

Detailed descriptions of equipment design and construction are 
available in Reference 3. Only a brief overview is presented here. 

As seen in Figure 1, the assembly consisted of a vertical test sur- 
face which contained a single column of eight rectangular protruding 
heated elements. Each protrusion simulated a 20-pin dual-in-line 
package. The elements were stainless steel blocks spaced one inch 
apart, center to center, on a plexiglass substrate. A 12.70 mm thick 
plexiglass shroud, when used, simulated the reverse side of a second 
printed circuit board. 

The blocks (Figures 2 and 3) had imbedded thermocouples on 
each of the five surfaces in contact with the fluid to provide surface 
temperature measurements. A sixth thermocouple was placed in the 
center of each element mounting slot for measurements of heater 
temperature. A miniature thermocouple probe mounted on a three- 
dimensional traverse enabled measurements of temperature distribu- 
tions within the fluid. 

A foil heater mounted on the component side facing the substrate 
was powered using a regulated D.C. power supply. Each heater was 
run in series with a precision resistor and all eight heaters were in 
parallel with the power supply. The current to each heater was 
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Figure 1 . System Configuration 
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0.02 R (4-PLACES ) 





Figure 2. Heater Block Schematic 
(after Hazard, Ref. 3) 
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12 ?6 20 2 ~ 26 



Figure 3. Mounted Foil Heater. Two Perspectives 

(after Hazard, Ref. 3) 
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calculated by subtracting the heater voltage from the source voltage, 
then dividing by the resistance of the precision resistor. Power dis- 
sipation through each element was determined by calculating the 
product of heater voltage and heater current. 

Figure 2 shows the orientation of the coordinate axes used in the 
discussion to follow. The front face of the component is that surface 
in the x-z plane in contact with the fluid. The side faces are the two 
x-y surfaces, right and left as determined by viewing the component 
from the front. The top and bottom faces are those in the y-z plane, 
downstream and upstream of buoyant fluid flow, respectively. 
Component height is in the x-direction, depth in y, and breadth in the 
z-direction. 

The immersion fluid was purified water contained in a one cubic 
meter, plate glass tank. Suspended in the water were particles of 
Pliolite, an inert pigment with a specific gravity of 0.93. The particles 
were illuminated with an eight milliwatt helium neon laser split into a 
plane by use of a cylindrical lens. The plane of visualization was varied 
by realignment of the laser lens assembly and the camera. 

B. EXPERIMENTAL PROCEDURE 

1. Established Method 

Details are again provided in Reference 3, hereafter referred 
to as the previous study. Noteworthy differences are noted here. 

2. Instrument Settings 

Power settings of 0.2, 0.5, 1.0, 1.2, and 1.5 watts were used 
for each channel width. Because of the inability to de-aerate the 
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immersion fluid, 1.5 watts was the upper limit of power available. 
Settings above this generated bubbles on the block surfaces. 

3. Instrument Readings 

When measuring steady-state surface temperatures, succes- 
sive temperature measurements at 10-minute intervals were com- 
pared. Only when all thermocouples varied less than 0.10° C were the 
final measurements taken. The temperature acquisition program of 
the previous study, modified to reflect the numbering of elements 
from bottom to top (flow direction) and to allow data storage 
(Appendix A), was then executed. 

When using the traversing thermocouple probe, the data 
acquisition system was programmed to continuously monitor the volt- 
age registered by a single channel. Observation of the on-screen print- 
out allowed determination of temperature variation from point to point 
within the flow field, as well as at any point over time. 

A program (Appendix A) which monitored from one to eight 
thermocouples over time was used when making a record of the tran- 
sient temperature response of the surface thermocouples to power 
step up or step down. Typical data presented were for the front face 
temperatures of the various components. 

4. Photographic Technique 

Best results were obtained when using an f stop of 2.8 and 
exposure time of 20 to 40 seconds. Camera placement was approxi- 
mately six inches from the tank wall at a height which allowed all 
eight blocks in the field of view. The camera was inclined at an angle 
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of about 5 degrees from the test surface when a shroud was present. 
This was helpful in minimizing reflected light. A shroud spacing of 
two times component height from the plexiglass substrate was deter- 
mined to be the minimum at which useful pictures could be obtained. 

C. DATA ANALYSIS 

While Qcond and Qconv were calculated exactly as recommended 
by the previous study, further data reduction differed significantly. 

For each block, a surface -averaged temperature was generated by 
multiplying the temperature and area of each fluid-exposed face and 
dividing the summation over the block by the total exposed face area 
of the block. The average over the eight blocks of the surface-averaged 
temperatures was then averaged with the immersion fluid tempera- 
ture to produce a film temperature. Fluid properties based on curve 
fits (Appendix B) of data from Incropera and Dewitt [Ref. 9] were then 
calculated for the film temperature. 

A new characteristic length, L , for determining a nondimen- 

Q 

sional temperature excess, -r-, a modified Grashof number, Gr*, and 

% 

Nusselt number, Nu, was chosen. The progression of its use for each 
block is as follows: 

5 

EXPOSED FACES 

£■ = - 4 

Y P 

^ EACH FACE 
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The various symbols used are defined in the List of Symbols section. 

Other characteristic lengths. L to determine block spacing along 
the plexiglass substrate and D to represent the horizontal breadth of 
the element, were also chosen. Their use is as follows: 
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-y- = nondimensional location of component centers along the 

test surface 

= horizontal position along the element’s front face 

Quantitative development was best expressed as plots of: 

• Nu vs. Gr* relationships 

• nondimensional block surface temperature excess vs. element 
position 

• temperature distributions in the fluid adjacent to the test surface 

• block surface temperatures vs. time 

Pertinent trends in these variations are discussed next. 
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HI. NO SHROUD- INFINITE CHANNEL WIDTH 



A. ALL ELEMENTS POWERED, STEADY STATE 
1. Flow Visualization 

Photographs of the natural convection flow in a plane (x-y) 
geometrically centered through each component are shown in Fig- 
ure 4. The visualizations represent all elements heated to power 
settings of 0.2, 0.5, and 1.0 watts, left to right, respectively. The 
exposure times were 40, 20, and 20 seconds, left to right. The pho- 
tograph toning techniques required for mass printing of this work do 
not allow sufficient resolution for detailed study. Should the reader 
desire to see prints of much greater clarity, contact Y. Joshi at the 
address given in the distribution list of this study. 

The flow appears to consist, in each case, of strong upward 
motion driven by the buoyant force near the test surface. The fluid 
tends to follow the contours of the protrusions. The particle traces 
are closer passing the element front face as compared to the inter- 
block space. Dead zones are not obvious except perhaps for very small 
areas where the blocks’ top faces intersects the substrate and for a 
slightly larger area above the top face of the last block in the array. 
There is no evidence of recirculation or formation of dead zones near 
the bottom face of any block. 
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Figure 4. Steady Flow in X-Y Plane at Power Levels of 0.2, 0.5, 
and 1.0 Watts. Powers Increase From Left to Right 
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In all three photographs, a developing boundary layer characteristic is 
evident farther out from the protrusions. It attains a thickness of 
about 3.5 component depths for 0.2 watts and 2 for 1.0 watt. Its 
origin moves upstream with increased power, from about 1 com- 
ponent height at 0.2 watts to 3 component heights at 1.0 watt. The 
velocities appear greater and the boundary layer becomes thinner as 
power is increased. Ar. ent fluid is entrained over the entire length 
of the layer. 

Figure 5 presents flow visualization in the plane (x-z) passing 
through the front faces of the elements. All elements were powered 
to 1.2 watts and exposure time was 40 seconds. The photograph 
shows an upward-moving buoyant layer near the components. Ambient 
fluid is entrained from each side into the flow, resulting in an increase 
in the horizontal extent of the upflow region down stream. The flow in 
the vicinity of the elements appears three-dimensional. 

As with the side view, the flow appears to follow the contours 
of the projections to about the same degree (dipping in about 3-5 mm 
from the side faces), the flow also appears to wrap around the cor- 
ners. The particle traces about 20 percent along the block (z direc- 
tion) from the side edge appear to be largely in the x-y plane and 
those at the side face are in the x-z plane. The particle traces in 
between are 3-D, each progressively moving along the 90° arc. 

Overall, the flow can best be described as following the con- 
tour of an hourglass split in the x-z plane. The elements are the upper 
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Figure 5. Flow in X-Z Plane. 1.2 Watts 



19 




and lower vessels. This configuration is repeated over the length of 
the array. It becomes more pronounced as the flow moves up the 
array and is eventually shed as a plume above the last element. 

2. Quantitative 

Figures 6 and 7 are nondimensional representations of tem- 
perature excess for 0.2 and 1.0 watts, measured on various component 
surfaces vs. block position. Also shown are the data for heater 
temperature measurements. The two plots follow the same patterns. 
The numbers for the 0.2 watts are, as expected, larger due to the 
inverse relationship with T*, which is h dependent. Also as antici- 
pated, the heater excesses are much greater. The drop in the heater 
temperature between blocks 6 and 7 could perhaps result from a 
slightly smaller (-4% less) power achieved by block 7. Unfortunately, 
the heater thermocouple for block 8 was inoperative. For both power 
settings, the imbedded surface thermocouples showed increasingly 
greater excesses along the array. 

Of primary importance in thermal design is the availability of 

the appropriate heat transfer relationships. As discussed earlier, such 

correlations for discrete protruding heat sources in liquids have not 

been available in literature up to now. To this end, a logarithmic plot 

» 

of Nu vs. Gr* (Figure 8) was generated with eight curves representing 
element position. Each curve had five points for the power settings 
0.2, 0.5, 1.0, 1.2, and 1.5 watts resulting in a range of Gr*. Although a 
linear pattern with Gr* was evident, the spread in data indicated an 
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Figure 6. Nondimensional Temperature Excess vs. 
Position. 0.2 Watts 
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Figure 7. Nondlmensional Temperature Excess vs. 
Position. 1.0 Watt 
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Figure 8. Data for Nusselt Number vs. 
Flux-Based Grashof Number 
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additional dependence on element location. A review of the data sug- 
gested a method by which it could be collapsed. A modified Nusselt 
number, Nu*, was calculated by use of a scaling factor based on the net 
convected energy of the flow at a downstream location. Since each 
block dissipates approximately an equal amount of power, the net 
convected energy is proportional to the number of upstream compo- 
nents. The scaling factor was of the form Bn a where Bn is the block 
number and a > 0. A value of a - 1/6 was found to provide the best 
correlation for the data. As expected, the strong three-dimensional 
effect on transport results in a weaker downstream dependence than 
the 1/5 exponent in the Fujii and Fujii correlation for a uniform flux, 
semi-infinite vertical surface. Therefore, 

l 

Nu*= Nu ■ Bn* 

for each element at each power setting. The modified data are plotted 
in Figure 9 along with the linear least squares fit. Based on this, the 
following correlating equation for the heat transfer rates is obtained: 

Nu*= L88(Gr *)' 5 

for 50 < Gr‘ < 1000 and Pr ~ 7. 
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Figure 9. Nusselt Number vs. Flux-Based Grashof Number 

With Least Squares Fit 
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B. ALL ELEMENTS POWERED, TRANSIENT 



1. Flow Visualization 

Figure 10 is a series of photographs taken during the 0.2 watt 
step-up transient. For each picture, the exposure time was 20 sec- 
onds. Exposure of the first picture began with power application. The 
others were initiated at 20, 80, 120, and 200 seconds, respectively, 
following the start of heating. Clearly, the first two show predom- 
inantly parallel flow. The next pictures represent increasing devel- 
opment of the contour flow and boundary layer characteristic 
development. By 200 seconds, almost fully developed flow is 
recognized. 

Figure 1 1 is a similar series, again of 20-second exposures, 
taken during the step-down transient from 1.0 watt. The first visual- 
ization at far left began at power down. The remaining photographs 
were initiated 80, 160, 240, and 440 seconds, respectively, following 
shutdown. The steady-state flow characteristics persist to an 
increasingly lesser degree for the entire period observed. The 
entrainment of ambient fluid lessens at large times. Additionally, flow 
origin moves slightly downstream over time. 

2. Quantitative 

The 1.0-watt step-up, step-down transient temperature 
responses for the component front surfaces are presented on the 
same graph (Figure 12). Steady state is achieved about 25 percent 
faster for the power up. The curves representing blocks 1 through 8 
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Figure 10. 0.2 Watt Step-Up Transient During 
Various Time Intervals 



27 











Figure 1 1 . 1.0 Watt Step-Down Transient During 
Various Time Intervals 
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Figure 12. Transient Temperature Response Following 
Initiation of Heating. 1.0 Watt 
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TIME CSEC> 



excluding 3 travel closely together, especially in the first quarter of 
the response period for the step. The almost single line corresponds 
to the parallel flow pattern seen in the photographs. This indicates a 
purely diffusive transport response for a large part of the transient. 
Later, convective effects become appreciable as entrainment develops 
downstream. For the step down, steady state temperature is reached 
much more quickly than breakdown of the established flow. The per- 
sistence of flow is conducive to the rapid cooling of the blocks. 

The plot for 0.2 watts (Figure 13) reaches steady state at 
about the same time for both step up and step down. This period is 
longer than for either 1.0-watt case and indicates that the stronger 
flow more readily establishes the final thermal pattern. 

C. TEMPERATURE MEASUREMENTS IN THE FLUID 

Characteristics of the thermal transport in the fluid were mea- 
sured using a traversing thermocouple probe adjacent to all protruding 
elements and for points midway between the components (Figures 14 
and 15). The fluid temperature responses near the blocks and inter- 
block spaces follow the same pattern, with temperature levels in the 
vicinity of blocks 4 and 6 surprisingly greater than for the others. 
Indeed, both graphs can be seen to have four groupings. Block 1 by 
itself, the lowest curve, is followed by 2, 3, and 5; then 4, 7, and 8; 
and finally block 6 by itself. More detailed local measurements are 
needed to explain these variations. 
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Figure 14. Fluid Temperature Distribution in the Normal 
Direction. Awav From Block Faces 
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Figure 15. Fluid Temperature Distribution in the Normal 
Direction. Awav From Plexiglass Substrate 
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The layer between the blocks is slightly thicker with respect to 
that at the block faces. This might be expected from the flow 
visualization (Figure 4) showing that the particle traces spread out 
between elements. In each case, the thermal layer is about one 
component height or less and is therefore deeply imbedded within 
the momentum exchange region. This is expected due to the 
relatively large Pr ~ 7 for water. 

Next, the traversing probe was placed in the thermal layer, 1 mm 
from the surface, and moved in the x-z plane across the test surface 
face (Figures 16 and 17). The temperature variations, with position, 
reinforce the photographs in suggesting strongly three-dimensional 
transport. Ambient temperature is the lowest point shown. The 
peaks at about 15 percent from each edge (Figure 16) correlate with 
the regions in the visualization (Figure 5), where the flow begins to 
shift from x-y dominant to x-z and may indicate a common region 
through which those particles are flowing. The result is more, 
warmer fluid convected over that region. Likewise, the peaks of Fig- 
ure 17 might represent a similar flow shift at locations between the 
components. 

Finally, the probe was placed above the center of the top face of 
block 8, 2 mm out from the substrate to measure the downstream 
temperature decay (Figure 18). The extremely high temperature in 
the immediate vicinity of the top face suggests a dead zone, as did the 
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Figure 16. Spanwise Fluid Temperature Distribution 
Across Components 



35 













l 



Figure 17. Spanwise Fluid Temperature Distribution 
Across Intercomponent Spaces 
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Figure 18. Fluid Temperature Decay Downstream 
of Top Component 
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visualizations, where a small pool of fluid is heated primarily by ther- 
mal conduction from the component upper surface. There was no 
sign of a vortex. 

All measurements with the probe were remarkably consistent 
over time (±1 pv) at any position. No evidence by any measurement 
indicated the presence of turbulent flow for the power levels 
investigated. 

D. INDIVIDUAL ELEMENTS POWERED 

Individual blocks were powered to 1.0 watt, one at a time. The 
resulting steady-state array temperatures were monitored and the flow 
visualized. These are discussed next. 

1. Flow Visualization 

A series of photographs show single-element heating as that 
element changes from one block to the next up the array. The expo- 
sure times for all pictures presented were 40 seconds. Figure 19 pre- 
sents left to right blocks 1, 3, and 4 heated and Figure 20 shows 
blocks 5, 6, and 7. The flow with only the lowest block heated resem- 
bles the all-components-heated case, indicating flow around block 
contours and a boundary layer characteristic over the entire length. As 
progressively higher blocks are heated, a quiescent area begins to 
develop adjacent to the lower part of the test surface. Typically, fluid 
entrainment results in vigorous flow around the first block upstream of 
the one heated. At two blocks upstream, a comparatively weaker flow 
is seen and at three the fluid is nearly quiescent. 
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Figure 19. Steady Flow Responses With Blocks 
1. 3. and 4 Individually Heated 
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Figure 20. Steady Flow Responses With Blocks 
5. 6. and 7 Individually Heated 
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2. Quantitative 



Nondimensionalized temperature excesses for the faces and 
heater of each element are displayed (Figures 21 through 25). The 
plots show very similar temperature excess patterns around the 
heated element as the selected component progresses up the array. 
The protruding elements have minimal effect when they are part of 
the unheated starting length. 

E. COMBINATIONS OF ELEMENTS HEATED 

As in Section D. progression was made up the array, except that 
all upstream blocks remained heated. Figures 26 through 30 indicate 
very little effect of powering of downstream components on the tem- 
perature responses of already-powered upstream components. Also 
evident is a very sharp decay in temperature in the unheated blocks. 
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Figure 21. Nondimensional Temperature Excess Levels at 
Front Faces of Individually Heated Components 
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Figure 22. Nondimensional Temperature Excess Levels at 
Right Faces of Individually Heated Components 
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Figure 23. Nondimensional Temperature Excess Levels at 
Top Faces of Individually Heated Components 
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Figure 24. Nondimensional Temperature Excess Levels at 
Bottom Faces of Individually Heated Components 
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Figure 25. Nondimensional Temperature Excess Levels at 
Heaters of Individually Heated Components 
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Figure 26. Nondimensional Temperature Excess Levels at 
Front Faces of Multiple Heated Components 



47 




Figure 27. Nondimensional Temperature Excess Levels at 
Right Faces of Multiple Heated Components 
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Figure 28. Nondimensional Temperature Excess Levels at 
Top Faces of Multiple Heated Components 
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Figure 29. Nondimensional Temperature Excess Levels at 
Bottom Faces of Multiple Heated Components 



50 




Figure 30. Nondimensional Temperature Excess Levels at 
Heaters of Multiple Heated Components 
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IV. SHROUD IN POSITION- FINITE CHANNEL WIDTHS 



A plexiglass plate (shroud) was used (Figure 1) to simulate the 
channel formed by two adjacent parallel printed circuit boards. The 
test surface as before has protruding heated elements. The shroud 
simulates the reverse side of the second board in an array. 

A. ALL ELEMENTS POWERED, STEADY STATE 
1. Flow Visualization 

Photographs presenting natural convection flow in the x-y 
plane, geometrically centered through each component (Figure 31), 
were evaluated. The visualizations are for all elements heated to 0.5, 
1.0, 1.2, and 1.5 watts. The exposure times were 20 seconds. The 
shroud distance was 18 mm from the substrate or approximately three 
component depths. The first six blocks are clearly visible in each pic- 
ture. Figure 32 shows the top portion of the test surface for 0.5 watts. 
The other power settings do not differ significantly and are not shown. 

The flow near the elements appears similar to that observed 
in the absence of the shroud. However, the shrouding wall prevents 
ambient fluid from entrainment into the buoyant up-flow in this plane. 
As power is increased, the paths described by the particle traces are 
more closely spaced, indicative of increased velocity and increased 
fluid flow. 
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Figure 31. Steady Flow in X-Y Plane. 0.5. 1.0. 1.2. and 1.5 Watts. 
18 mm Channel Spacing. Bottom Blocks 
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Figure 32. Steady Flow in X-Y Plane. 0.5 Watts. 18 mm 
Channel Spacing. Tod Blocks 
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Figure 33 is the front view of flow in the x-z plane through 
the front surface of the elements. It is for 1.0 watt, 18 mm shroud 
spacing, and exposure time of 40 seconds. Unlike the no-shroud con- 
dition, flow is more restricted to the sides of the protruding blocks. 
The characteristic hourglass flow seen for infinite channel width is 
weak in the vicinity of the middle two blocks and absent elsewhere. 
The entrained fluid turns very sharply into the flow direction. Because 
the flow is restricted to the channel width and is predominant on the 
block sides, its velocity is much greater here. Consequently, there is 
more forceful entrainment of fluid in the x-z plane. 

2. Quantitative 

Due to the large amount of data collected, temperature 
excesses for power settings of 1.2 and 0.2 watts only are presented in 
graphical form. These settings are representative of the upper and 
lower ends of the power range available. Channel widths presented 
are 6.2, 8, 14, 18, and 23 mm measured from the substrate. 

The first width was just barely greater than component 
height. This spacing is representative of pure conduction dominated 
transport. The nondimensionalized temperature excess vs. element 
position graphs for 0.2 watts (Figures 34-38) and 1.2 watts (Figures 
39-43) show significant difference between the smallest spacing and 
all others. With that spacing excluded, Figures 44 through 53 present 
the relevant data in a manner better suited for analysis. The plots for 
0.2 watts show the next closest spacing, 8 mm from the substrate, to 
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Figure 33. Flow in X-Z Plane. 1,0 Watt. 18 mm Channel Spacing 
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Figure 34. Nondimensional Temperature Excess. 0.2 Watts. 
Various Spacings. Front Faces 
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Figure 35. Nondimensional Temperature Excess. 0.2 Watts. 
Various Spacings. Right Faces 
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Figure 36. Nondimensional Temperature Excess. 0.2 Watts. 
Various Spacings. Top Faces 
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Figure 37. Nondimensional Temperature Excess. 0.2 Watts. 
Various Spacings. Bottom Faces 
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Figure 38. Nondimensional Temperature Excess. 0.2 Watts. 

Various Spacings. Heaters 
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Figure 39. Nondimensional Temperature Excess. 1.2 Watts. 
Various Spacings. Front Faces 



62 




Figure 40. Nondimensional Temperature Excess. 1.2 Watts. 
Various Spacings, Right Faces 
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Figure 41. Nondimensional Temperature Excess. 1.2 Watts. 
Various Spacings. Tod Faces 
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Figure 42. Nondimensional Temperature Excess. 1.2 Watts. 
Various Spacings. Bottom Faces 
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Figure 43. Nondimensional Temperature Excess. 1.2 Watts. 

Various Spacings. Heaters 



66 




Figure 44. Nondimensional Temperature Excess. 0.2 Watts. 

Excl. 6.2 mm. Front Faces 
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Figure 45. Nondimensional Temperature Excess. 0.2 Watts. 

Excl. 6.2 mm. Right Faces 
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Figure 46. Nondimensional Temperature Excess. 0.2 Watts. 

Excl. 6.2 mm. Top Faces 
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Figure 47. Nondimensional Temperature Excess. 0.2 Watts. 
Excl. 6.2 mm. Bottom Faces 
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Figure 48. Nondimensional Temperature Excess. 0.2 Watts. 

Excl. 6.2 mm. Heaters 
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Figure 49. Nondimensional Temperature Excess. 1.2 Watts. 

Excl. 6.2 mm. Front Faces 
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Figure 50. Nondimensional Temperature Excess. 1.2 Watts. 

Excl. 6.2 mm. Right Faces 
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Figure 51. Nondimensional Temperature Excess. 1.2 Watts. 

Excl. 6.2 mm. Top Faces 
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Figure 52. Nondimensional Temperature Excess. 1.2 Watts. 
Excl. 6.2 mm. Bottom Faces 
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Figure 53. Nondimensional Temperature Excess. 1.2 Watts. 

Excl. 6.2 mm. Heaters 
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yield higher temperatures over the range of element positions. For 
higher spacings, there was no consistency in the relative position of 
the curves from one element position to the next. This indicates 
relative independence of transport on channel width for spacings 
greater than 8 mm. However, data for 1.2 watts includes even the 
8 mm spacing within the mix of curves (its curve does not consistently 
yield higher temperatures). Therefore, for 1.2 watts the indepen- 
dence of transport extends below the 8 mm spacing. General patterns 
observed for both powers are: 

• a steady rise in front and side face temperatures with element 
position as expected due to exposure to increasing convected 
energy at higher positions. 

• a relative levelling of top and bottom temperatures at higher posi- 
tions. This is perhaps due to increasing fluid velocities along the 
channel. 

• the drop in heater temperatures at block 7, which has four per- 
cent less power on average. 

The differences between the two powers result from the 
stronger driving force at 1.2 watts. Increased fluid velocities result in 
improved heat transfer characteristics, making 8 mm an acceptable 
spacing. This point is emphasized by referring to the graphs of com- 
ponent face temperature excess vs. shroud spacing for 0.2 watts 
(Figures 54-58) and 1.2 watts (Figures 59-63). Significantly, for the 
1.2 watts there is no further temperature reduction for a spacing 
greater than 8 mm. Indeed, there is a slight rise in temperature of 
the latter blocks (7 and 8) when the shroud is removed. These higher 
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Figure 54. Nondimensional Temperature Excess vs. Shroud 
Position. 0.2 Watts. Front Faces 
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Figure 55. Nondimensional Temperature Excess vs. Shroud 
Position, 0.2 Watts. Right Faces 
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Figure 56. Nondimensional Temperature Excess vs. Shroud 
Position. 0.2 Watts. Top Faces 
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Figure 57. Nondimensional Temperature Excess vs. Shroud 
Position. 0.2 Watts. Bottom Faces 
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Figure 58. Nondimensional Temperature Excess vs. Shroud 
Position. 0.2 Watts. Heaters 
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Figure 59. Nondimensional Temperature Excess vs. Shroud 
Position. 1.2 Watts. Front Faces 
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Figure 60. Nondimensional Temperature Excess vs. Shroud 
Position. 1.2 Watts. Right Faces 
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Figure 61 . Nondimensional Temperature Excess vs. Shroud 
Position, 1.2 Watts. Top Faces 
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Figure 62. Nondimensional Temperature Excess vs. Shroud 
Position. 1.2 Watts. Bottom Faces 
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Figure 63. Nondimensional Temperature Excess vs. Shroud 
Position. 1.2 Watts. Heaters 
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temperatures (at the last blocks) are not observed when a shroud is in 
position because the global forced convection within the channel is 
driving the fluid more strongly past the last block. As seen from Fig- 
ures 54 through 58, face temperatures at a power level of 0.2 watts for 
8 mm are somewhat higher than for greater spacings. By 14 mm, the 
difference with larger spacings begins to become indistinguishable. 
Clearly for this geometry and for the power settings examined, a 
channel width of 2.5 times component depth is most satisfactory 
while a width of 1.3 times element depth may be acceptable. Further 
confirmation of the observation that spacings beyond 8 mm do not sig- 
nificantly affect heat transfer is seen in Figures 64 through 68. These 
show the Nu* vs. Gr* data along with the correlation, Nu* = 1.88 
(Gr*) 15 , obtained in Chapter III. The correlation for an unshrouded 
surface holds very well through 14 mm. At 8 mm there is a slightly 
greater deviation. As anticipated, it is inapplicable for the smallest 
spacing of 6.2 mm. 

B. INDIVIDUAL ELEMENTS POWERED 

As in the absence of a shroud, individual blocks were powered to 
1.0 watt, one at a time, but with a shroud spaced 14 mm from the sub- 
strate. The resultant data is displayed in Figures 69 through 73. The 
plots follow similar trends as for the no-shroud condition (Figures 21- 
25) but show somewhat smaller temperatures. Again, the protruding 
elements have minimal effect on the unheated starting length. 
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Figure 64. M/* vs. Gr*. With Correlation. 23 mm Channel Spacing 
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Figure 65. Nu* vs. Gr*. With Correlation. 18 mm Channel Spacing 
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Figure 66. Nj£ vs. Gr*. With Correlation. 14 mm Channel Spacing 
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Figure 67. Ni£ vs. Gr*. With Correlation. 8 mm Channel Spacing 
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Figure 68. Nul vs, Gr*. With Correlation. 6.2 mm Channel Spacing 
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Figure 69. Nondimensional Temperature Excess. Front Faces. 
Individually Heated Components. 14 mm Channel Spacing 
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Figure 70. Nondimensional Temperature Excess. Right Faces. 
Individually Heated Components. 14 mm Channel Spacing 
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Figure 71. Nondimensional Temperature Excess. Top Faces. 
Individually Heated Components. 14 mm Channel Spacing 
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Figure 72. Nondimensional Temperature Excess. Bottom Faces. 
Individually Heated Components. 14 mm Channel Spacing 
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Figure 73. Nondimensional Temperature Excess. Heaters. 
Individually Heated Components. 14 mm Channel Spacing 
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C. COMBINATIONS OF ELEMENTS HEATED 

In a manner identical to that without a shrouding wall, progres- 
sion of elements heated was made up the array while all upstream 
blocks remained powered. Additionally, block 8 was powered 
throughout and a shroud 14 mm from the substrate was in place. The 
plots (Figures 74-78) reveal that for this channel width, as with no 
shroud, the powering of downstream components has minimal effect 
on those upstream components already powered. The plots further 
show a weak effect by a heated element (block 8) on those compo- 
nents between it and the ones powered upstream, which are heated 
only by the convective flow generated upstream of them. 

D. ALL ELEMENTS POWERED, TRANSIENT 
1. Flow Visualization 

Figure 79 shows transient flow development in the x-y plane 
following 1.0 watt power step-up with a 19 mm spaced shroud in 
place. All exposures are for 20 seconds. The photographs from left to 
right are initiated at 0, 20, 40, 100, and 140 seconds following start of 
heating. In the first picture, parallel flow is developing over blocks 3 
through 8 with relatively little movement at blocks 1 and 2. By picture 
number 3, the flow follows more closely the contours of the protrud- 
ing blocks and originates further upstream. By photograph 5, the flow 
appears nearly fully developed. 

The power step-down from 0.2 watts was again for a 19 mm 
channel width and 20-second exposures (Figure 80). The pictures 
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Figure 74. Nondimensional Temperature Excess. Front Faces. 
Multiple Components Heated. 14 mm Channel Spacing 
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Figure 75. Nondimensional Temperature Excess. Right Faces. 
Multiple Components Heated. 14 mm Channel Spacing 
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Figure 76. Nondimensional Temperature Excess. Top Faces. 
Multiple Components Heated. 14 mm Channel Spacing 
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Figure 77. Nondimensional Temperature Excess. Bottom Faces. 
Multiple Components Heated. 14 mm Channel Spacing 
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Figure 78. Nondimensional Temperature Excess. Heaters. Multiple 
Components Heated. 14 mm Channel Spacing 
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Figure 79. Flow Response During 1.0 Watt Step-Up 
Transient. 19 mm Channel 
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Figure 80. Flow Response During 0.2 Watts Step-Down 
Transient. 19 mm Channel Spacing 
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from left to right were initiated at 0, 60, 200, 320, and 460 seconds, 
respectively, following cessation of power. As with infinite channel 
width, the flow characteristic persists over the entire period observed. 
Only in the last picture does the flow begin to follow less closely the 
contour of blocks 7 and 8 and originate farther downstream. 

2. Quantitative 

The transient temperature response curves for powers of 1.0 
and 0.2 watts and a shroud spacing of 19 mm (Figures 81 and 82) fol- 
low exactly the same pattern as those for no shroud (Figures 12 and 
13). For the spacing examined, the shroud does not affect the trans- 
port evolution. 
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Figure 81. Transient Temperature Response 
Following Initiation of Heating. 1.0 Watt 
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Figure 82. Transient Temperature Response 
Following Cessation of Heating. 0.2 Watts 
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V. CONCLUSIONS 



Flow visualization in concert with quantitative analysis provided a 
number of significant determinations: 

• The convective flow was three-dimensional. 

• The flow was laminar. 

• The thermal layer was imbedded within the momentum exchange 
region. 

• Heated elements have minimal effect on those elements, 
upstream, which are already heated. 

• There is no gain in heat transfer for channel widths greater than 
two times component depth. 

• A correlation has been found to predict heat transfer rates for the 
subject geometry. 
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VI. RECOMMENDATIONS 



It is suggested that the following areas of study be experimentally 
explored: 

• Using the same test surface: 

- explore transient response for other block faces 

- with the traversing thermocouple probe, determine variation of 
temperature with position (in the thermal layer) for all blocks 

• Construct new test surface and rerun the same battery of experi- 
ments to determine significant changes when: 

- heated elements are flush mounted 

- heated elements are spaced closer in the x-direction. 



Ill 



APPENDIX A 



SOFTWARE 

TEMPERATURE ACQUISITION PROGRAM 



10 n i l u H i i M ! ! ! i M M M ! ! i M 1! I M ! I M ! ! !! ( M M ! * M ! I M ! i • 

20 I' TEMPERATURE ACQUISITION PROGRAM 

30 ' ! ( STEAOV STATE ) 

40 M M l M M i M M i I | I I | i M M M i I I I l j M i M ! I M j ! I I I I * ! ! i ! I I I 



50 


REAL Vol t s< 60 ) 




60 


REAL Temp< 59 ) 




70 


CREATE BOAT ” THES 


: ,700,0,2 M .48 


80 


ASSIGN 6Path3 TO "THES : ,700,0,2" 


90 


PRINT " 


BLOCK #8'* 


100 


PRINT M " 




1 10 


OUTPUT 709 ; "CONFMEAS 


DCV, 100, USE 0” 


120 


ENTER 709 ; Vol t s( 60 ) 




1 30 


OUTPUT 709 ; '‘CONFMEAS 


DCV, 100- 105. USE 0 M 


140 


FOR 1-0 TO 5 




150 


ENTER 709;Volts( I ) 




160 


Tenp( I >-.0006797+< 25825. 1 328»Volta( I ) )-( G07789 . 24G7»( Vo lta(I)» Volt a< !)))-( 



21 352034. 3364*< Vol ts< I ) * 3) > + < 83708 10996. 1 874* ( Volt s< I ) A 4) ) 



170 


PRINT “T. C. # M j 1+1 .“ 


Volts 


D.C. " ; Vol ts( I ) , " Tenp. DEG. C M ;Temp(I) 


180 


NEXT I 






1 90 


PRINT " ” 






200 


PRINT ” 




BLOCK 47" 


210 


PRINT " “ 






2 20 


OUTPUT 709;" CONFMEAS DCV 


. 106- 


1 l 1 .USE 0" 


2 30 


FOR 1-6 TO 11 






240 


ENTER 709;Volts< I ) 






2 50 


T emp< I )-.0006797+( 2 5825. 


1 328* Vol t a( I ) )-( 607789. 24G7»( Vol ts< I ) • Vol tsU ) ) ) 



2 1952034. 3364* < Volt s< I r 3) ) + < 8370810996. 1874*< Vol ts< I ) "4 ) ) 

260 PRINT "T. C. ril + l." Volts O.C. ’* ; Vol t s( I ) , "Temp. DEG. C " ; Tenp< I ) 

270 NEXT I 

280 PRINT M M 

290 PRINT '• BLOCK 46" 

300 PRINT " " 

310 OUTPUT 709; "CONFMEAS OCV ,112-1 17, USE 0" 

320 FOR 1*12 TO 17 
330 ENTER 709;Volts( I ) 

340 Tewp< I ) * . 0006797+< 25825. 1 32 8* Vol t s< I ) ) -< 607789 . 2467* ( Vol ts( I ) * Vol ts< I ) ) >-< 
2 1952034. 3364* (Volt s< I > * 3 ) ) +< 8370810996. 1 874* ( Vol t s< I ) “ 4 ) ) 

350 PRINT "T. C. rihl," Volts O.C. " ? Vol ts< I > ."Temp. OEG. C " ; Temp< I > 

360 NEXT I 

370 PRINT M M 

380 PRINT M BLOCK 45" 

390 PRINT “ " 

400 OUTPUT 709 ; " CONFMEAS DCV , 1 1 8- 1 1 9 . USE 0 M 
410 FOR 1=18 TO 19 
420 ENTER 709; Volts< I ) 



l i i 
I i 
i i 
i i i 



4 30 Tenp< I > - . 00®6797+< 258Z5. 1 3Z8»Volts< I ) >-(607789 . 2467* < Volt s( I ) • Vol ts< I > > >-< 
2 I 952034. 3364* < Vol ts( I >“ 3) > + (8370810996. 1 874* ( Volt s< I > “4 ) ) 

440 PRINT "T. C. #";I+I," Volts O.C. “ i Vol ts( I > . "Temp. DEG. C " ; Tenp( I > 

450 NEXT I 

460. OUTPUT 709 i "CONFMEAS DC V , 200-203 . USE 0" 

470 FOR 1-20 TO Z3 
480 ENTER 709; Volt s( I > 

490 Tenp( I >-.0006797+1 Z58Z5. !3Z8*Volts( I > >-( 607789 . Z467*< Vol ts( I ) *Vol ts( I ) ) >-< 
21952034. 3364* ( Volt s< I >*3> > + ( 8370810996. 1 874* ( Volt s( I >*4) ) 



500 


PRINT M T. C. 


4"?I+1," Volts D.C. " ;Volts< I ) /’Tenp. DEG. C "?Temp<I> 


510 


NEXT I 




520 


PRINT M M 




530 


PRINT M 


BLOCK 44" 


540 


PRINT ” M 




550 


OUTPUT 709 ?’ 


•CONFMEAS DCV ,204-209 , USE 0" 


560 


FOR 1*24 TO 


29 


570 


ENTER 709 ; Vol t s( I ) 


580 


Tempi I >«. 0006797+1 25825. 1 328* Volts! I > >-( 607789 . 2467* ( Volt s< I > » Vol ts( I > > > 



2 1952034. 3364* ( Volts< I >* 3 >>+< 83708 1 0996 . I 874* ( Vol ts( I >*4>) 



590 


PRINT M T. C. 4"?I+1 


( M Volts D.C. M ;Volts( I ) , M Temp. DEG. C " ; Tenp< I ) 


600 


NEXT I 




610 


PRINT " M 




620 


PRINT M 


BLOCK 43" 


630 


PRINT M M 




640 


OUTPUT 709 ; n CONFMEAS 


DCV ,210-215, USE 0" 


650 


FOR 1*30 TO 35 




660 


ENTER 709 ? Vol t s< I ) 




670 


T emp( I )* . 0006797+ < 25825. !328*Volts< I ) )-< 607789 . 2467* < Vol t s< I ) * Vol t s( I ) ) ) 



2 1952034. 3364* ( Volt s( I >*3) > + < 8370810996. 1 874* ( Volt s( I 1*4) > 



680 


PRINT "T. C. 4" tl + 1 , 


, " Volts D.C. " ? Volts< I ) ."Temp. DEG. C M ;Tepip<I) 


690 


NEXT I 




700 


FOR J-l TO 14 


fc. 


710 


PRINT " " 


- 


720 


NEXT J 




730 


PRINT " 


BLOCK #2" 


740 


PRINT *' " 




750 


OUTPUT 709? "CONFMEAS 


DCV. 2 16-2 19. USE 0" 


760 


FOR 1*36 TO 39 




770 


ENTER 709 ; Vol t s( I > 




780 


Temp( I )-.0006797+( 25825. 1 328»Volts( I ) )-< 607789 . 2467* ( Vol t s( I >*Volts( !>>)-( 



21952034. 3364* (Volt s( I >*3) )+< 8370810996. 1874*< Vol ts< I > “4 ) > 

790 PRINT "T. C. #";I + 1." Volts D.C. " ; Vol ts( I > . "Temp. DEG. C " ; Tempi I ) 

800 NEXT I 

810 OUTPUT 709 ; "CONFMEAS DCV , 300- 301 . USE 0" 

820 FOR 1-40 TO 41 

830 ENTER 709;Volts(I> 

840 T enp( I > - . 0006797+( Z5825 . 132 8*V 0 lts( I > )-< 607789. Z467*( Volt s< I >*Volts< I > > >-< 

2 1952034. 3364* ( Vol ts( I >*3> > + (8370810996. 1874*< Volts< I >*4> > 

850 PRINT "T. C. JTiI+l." Volts D.C. " ; Vol ts< I > , "Temp. DEG. C ";Temp(I> 

860 NEXT I 

870 PRINT " " 

880 PRINT " BLOCK #1" 

890 PRINT " " 

900 OUTPUT 709 ; "CONFNEflS DCV . 302-307 , USE 0“ 

910 FOR 1-42 TO 47 

920 ENTER 709;Volts( I > 



1 13 



9 30 T enp( I )=. 0006797+1 25825. I 328+ Volts! I ) )-( 607789 . 2467* ( Vol t s( I > • Vol tsl I ) ) )-< 
2 1 952034. 3364M Volts! I ) A 3 ) > + < 83708 1 0996 . I 874»< Vol ts! I > A 4 > > 

940 PRINT "T. C. #";I+1." Volts D.C. " ; Vol t s< I ) . "Temp . DEG. C " i Temp! I ) 

950 NEXT I 
950 ! 

970 i 

980 OUTPUT ePath3; Temp! » > 

990 ! 

1000 t 

1010 PRINT " " 

1020 PRINT " " 

1030 PRINT " BATH TEMPERATURES < TOP TO BOTTOM)" 

1040 PRINT " " 

1050 OUTPUT 709 ; “CONFME AS DCV . 317-319 .USE 0" 

1060 FOR I >57 TO 59 
1070 ENTER 709! Volts( I ) 

1080 Temp! I >-.0006797+! 2 5825. 1 32 8* Vo 1 1 s ( I ) ) - < 607789 . 2467»( Vol ts< I >*VoIts< I ) > >-< 
2 1952034. 3364»( Volt s< I ) * 3) > + ( 8370810995. 1 874* < Volts< I > A 4> > 

1090 PRINT "VOLTS D.C. " ; Volts! I > ."TEMP. OEG. C "iTemp!I> 

1 100 NEXT I 
1110 ENO 
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10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

1 10 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

2 30 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 



DATA REDUCTION PROGRAM 



I ; i i ! I { ! M |l | ! I m j i i { 1 ! i U M 1 ! < ! M M ! ( M ! I M ! ! ! ! l ! ! I l M i ! i M l ! I M 1 

l i l ! 

M DATA REDUCING PROGRAM (BASIC) !! 

M ii 

| I ! i l i i ! i l M | M I M | i I M i M ! M ! ! M ! ! ! ! ! ! ! M M ! M ! ! ! ! M ! M M ! ! ! ! I i ! ! 

REAL Tenp( 48) , Tt< 48) ,Qcond< 8) f Qconv< 8) . Qf 1 ux ( 8 ) , H< 8 ) , Nu< 8 ) ,Nustar< 8) 

REAL Gr< 8) ,Ttt< 8) ,At< 8) , Dt f ( 8 ) , Dt r< 8 ) , Ot 1 < 8 ) ,Dtt< 8) ,Otb< 8) ,Dth( 8) ,B1< 8) 
REAL Jstar( 8 ) , Jp( 8 ) ,Grstar( 8) 

REAL Power , Space , Dxpg , Dxr , Kpg , Kr .Block 
REAL A.G.Dis , Try , Kh2o .Beta, Kinvisc.Pr 
INTEGER I ,J,K .L.M.N.O 
! 

Dx pQ* . 006731 
Dxr* . 003 1 75 
K pg- .1421 
Kr» . 0389 
Chpht-. 00787 
Chpdpth- . 006 1 0 
Chpudth- . 02 388 

Pnn«( 2*( Chpht +Chpudth ) )-*•< 2*( 2*< Chpht+Chpdpth ) ) ) + ( 2*( 2*( Chpudt h+Chpdpth ) ) ) 

A 1 -Chpudth* Chpht 

A2 -Chpudth* Chpdpth 

A3*Chpht * Chpdpth 

Atotal-A1+(2.*A2)+< 2.*A3> 

At otal a-( 2 . * A2 ) + ( 2 . ♦ A3 ) 

G=9 . 81 

I i i i i I | | I ; ! I i ! M ! I i i ! ! ! i I { M l j i ! !! | i j | I ! I ! | I I | i i | l 1 | ! I | j i I ! I | i l I I | l| i M l 

i! CALCULATE CHARACTERISTIC LENGTH M 

i I I l I I i I I | I I l i j! M l l I I i l i 1 I I I i I | I hi j! I I l i I I I i | i | I I 1 I i l I | i i I i i I i I i i ! 1 !! ii 
L bar-Atotal /Pr i n 

i 

R1 =Dxpg/( Kpg*A1 ) 

R2*Dxr/< Kr* A 1 ) 

i 

PRINTER IS 1 

PRINT "INPUT RUN tt . INTEGER ONLY" 

INPUT Try 
PRINT " " 

PRINT M M 

PRINT "INPUT POWER VALUE IN UATTS" 

INPUT Pouer 
PRINT " " 

PRINT " " 

PRINT "INPUT BLOCK S POUERED" 

INPUT BlnpS 
PRINT " " 

PRINT " " 

PRINT "INPUT SPACING IN MM" 

INPUT Space$ 

PRINT " ” 



510 PRINT " " 

5Z0 PRINT "INPUT AMBIENT TEMPERATURE" 

530 INPUT Tinf 

540 PRINT " •* 

550 PRINT " " 

560 l PRINT "INPUT NEU FILENAMES" 

570 ! INPUT NeufileS 

580 ! PRINT " " 

590 ! PRINT " " 

600 PRINT "INPUT OLD FILENAMES" 

610 INPUT 01 df i 1 e$ 

620 i 

630 ! 

640 ! CREATE BDAT NeufileS, 128 

650 i ASSIGN ©Road TO NeufileS 

660 ! 

670 ASSIGN ©Street TO OldfileS 

680 ENTER ©Street ; Tenp( * ) 

690 FOR 1-0 TO 47 

700 J-J+1 

710 IPRINT J , T enp( I ) 

720 NEXT I 

730 i 

740 ! 

750 ! i | I l I M I |i l M i I I M I I I I I M | I i | I M M M M M M • M ! !! ! M ! ! ! ! M M ! ! ! ! ! 

760 ! i CALCULATE AVERAGE SURFACE TEMPERATURES i i 

770 ! M ! M ! ! I ! I ! I! ! ! M ! M l i l ! M i M ! ! I ! I ! ! I! I ! ! ! M ! ! ! ! ! ! ! ! ! ! ! i ! I ! 1 ! ! ! ! 

780 ! PRINT "AVG. SURFACE TEMPERATURE- 

790 J-0 

800 FOR 1-1 TO 8 

810 IF 1-6 THEN 

820 At( I )*( ( Tenp( J + l )*A3)+( Temp(J + 2 )*A3>+( Tenp( J+3)*A2)+(Temp(J+4)*AZ))/ 

Atotala 

830 ELSE 

840 At ( I ) -( ( T enp< J ) * A 1 ) + < Temp< J+l ) * A3 > + « Tenp< J + Z >*A3)+< Tenp( J+3)*AZ ) + < Te 

np( J+4 ) * A2 ) )/Atotal 

850 END IF 

860 i PRINT At ( I ) 

870 J-J+6 

880 NEXT I 

890 i 

900 M i i i I i j i i I i i I i i i i I i M I l M l i i l i l i j i i i i i I | i i I i i i i i i i M i { i i i i i M i l l 

910 M CALCULATE TFILM i ! 

g 20 i i i i i i i i t i i i i ; i i i i i i i i i i i i I u i i i (I ; i ; i i i i i i i i i i i i i i i i i i i i i | i ii i i i 

930 J-l 

940 Tbar-( At( J)+At( J+l ) + At( J + 2 )+At( J+3)+At( J + 4)+At( J + 5 ) + At ( J + 6 ) +At ( J+7) )/8.0 
950 Tf ii nc=( ( Tbar+T i nf )/Z . 0 ) 

960 Tf llm-Tf ilmc+273. 15 

970 ! PRINT "TFILM DEGREES C",Tfilmc 

980 ! PRINT "TFILM" . Tf l In 

990 M l ( i ! M ! ! M i M i m i i i i I I i ii m m | M I M M | i M I l M I i I M i M m I i i l l l M 

1000 H CALCULATE UATER PROPERTIES U 



1010 

10Z0 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1 100 

1 1 10 

1 120 

1 130 

1 140 

1150 

1 160 

1170 

1 180 

1190 

1Z00 

1 Z 1 0 

1ZZ0 

1 Z 30 

1 Z40 

1 Z50 

1 Z60 

1270 

1 Z80 

1 Z90 

1300 

1310 

1 3Z0 

1 330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

14Z0 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 



! i ! I i i M M I I | I I l l | I M ; I ! I l | j I I l M M l l ! II I M I ! ! I t ! ! I l l | I i | | I ! i j i t { 

Bet a-( -7.9448548E-8M Tf 1 1 2 ) ) + < 5 . 7356479E- 5* Tf 1 1 m ) - . 00978 1 0563 
Spvol-<4.S9699E-9*<Tfilm A Z> )♦< -Z . 53745E-6*Tf il n )+ . 001 341903 
Dynvisc-< 3. Z 34851 1E-7*< Tf i 1 n*Z )> + <-. 000Z 1 474487*Tf 1 1 pi )+. 0361 6679Z 
K invisc-Spvol *Dynvisc 
KhZo-< 1 .4181 8181 8E- 3* Tfilm>+. 1866 

Pr-< 4 . 65706E- 3* ( Tf iln*2 ) >♦< -Z . 922094* Tf il* >+463. 3319 

il | { l ! I ! M il ! ! ! | ! I i I M ! ! ! ! ! ! ! 1 ! ! ! M ! ! ! M M ! ! ! ! ! ! M M ! M ! I ! ! i I! ! ! 1 

!! CALCULATE TEMP- TINF !! 

! I !!!! i !! M !!!!!!!!! I 1 M )!! I !!!!!! M !!!!! II I !!!!!! M !!!!!!!!! 1 I ! ! 
i PRINT ” TEMP- TINF" 

FOR 1-0 TO 47 
IF 1-5 THEN 
Tt( I >-888.888 
ELSE 

IF 1-30 THEN 
Tt < I >-888.888 
ELSE 

Tt< I )*Temp< I )-Tinf 
END IF 
END IF 

! PRINT Tt < I > 

NEXT I 

!!!!!! i !!!!!! 1! !!!! i !!! M !!!!! I 1! I !!!!!!!!!!!! I !!!!!! M !!! I !!!!!! ! 

! ! CALCULATE CONDUCTION LOSSES THROUGH THE TEST SURFACE ! ' 

I !!!! M !!! || !! i !!!!!!!!! I |M I | M I l !! Il !!!!!!!!!! I !!!!!!!!!!!!!! M ! 

! PRINT "QCOND" 

L-5 

FOR 1-1 TO 8 
IF 1-1 THEN 

Qcond( I ) • . 0 1 5*Pouer 
ELSE 

QconcK I >-Tt< L >•< 1 .0/< R1+RZ > > 

END IF 

i PRINT Ocond(I) 

L-L + 6 
NEXT I 

I ! it l | I I | | I i I I I ; i ! | | ; | i | | i i i | | i i ! | I i i i i ! | i I I I i | I | I I | | | M | i i | i i i I I I i 

» ! CALCULATE CONNECTED HEAT FLUX * ' 

i I | I | i I I ! i i I I ! I ! | i I I | | i I I | ! I I ! i i | I | i | I I | | | | i I i | | i I I I i i | i | | I I i I M | | i 

» PRINT " OCONV" 

FOR 1-1 TO 8 

Oconv( I )-Pouer-Qcond< I > 

Qf lux< I ) -Qconv< I > / A 1 
NEXT I 

I I I I | | I I I I | | I | | I I I | | I | I I | | I I | I I | I I I 1 I I I I | I I | I I M I I | I I M t I I I 1 I It M I I 

M CALCULATE TEMPERATURE SCALING FACTOR M 

I i | I | l I ! ! I ! I i I | I I l M ! ! I I M ! ! M ! I I ! M I I I I M M I I ! I I ! >! I I I I i I I I I I ! t I ! I 

! PRINT " T * " 

FOR 1-1 TO 8 
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188® 

189® 
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191® 

192® 
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195® 
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198® 
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Tstar< I ) -Qf lux ( I )*Lbar/Kh2o 
! PRINT TSTAR< I ) 

NEXT I 

i l M I i{ i i M l I l | I M ! M M I ! M ! I ; I i I I i i ! I I I I | I i ! M I M I I ! I < I I I M ! I I M t j 

! ! CALCULATE MODIFIED GRASHOF NUMBER ! ! 

i | l M i ! I i ! I M ll j M ll i | l M i I l I i l l l i i i ! i I i ! M i M ! M ! < M ! I < i i m I | I j | M 

FOR 1*1 TO 8 

Grstar( I )*G*Beta*Qf lux( I ) * ( Lbar A 4)/< Kh2o*< K invisc A 2 ) ) 

NEXT I 

i M i H M t ! M 1 M I i ! M i M * ! M < M < ! ! I I i ! H i M H !! ! M M M M ! i I i i l ! M M < 

M CALCULATE HEAT TRANSFER COEFFICIENT M 

j ! ! I ! i ( I ! M 1 l M M 1 M M !! i M H ! ! ! M i i H ' ! ! ! H !! M M !! ! M ! i ( < i M H ! ! ! 

FOR 1*1 TO 8 

H( I )*Qf lux< I )/< At ( I )-T inf ) 

NEXT I 

! M | M I ! M M 1 |l M I ! ! 1 ! ! ! ! I ! I M l l n i I l l j M i j | { i ! M ji j 1 | ji I | i i i M M l 

! CALCULATE NUSSELT NUMBER ! 

| M M i M M | M i i | I ! M | j ij M H i I I I M i M i ! I M ! M ! M I I i ! M I ! i ! 1 M ! ! ! ! ! 
FOR 1-1 TO 8 
Nu< I ) *H( I )*Lbar/Kh2o 
NEXT I 

l M I i II M ! ! ! ! I ■ I I M 1 ! I I ! I ! ! ! ! I I I M I ! ! I I I M I ! ! H I 1 I I I ! ! M I ! I M ! M M I 

1 CALCULATE NUSSELT SCALING FACTOR !! 

I i l l ! ll M i | I M M ! II m i m II I M I ll i I i M I I ! 1 ! ! ! ! I ! i ! I M ! ! ! I I I I ! I I M I M 

FOR 1-1 TO 8 
Js*9 - 1 

Nust ar ( I )-Nu< I )*Js A < 1/6) 

NEXT I 

! I i | | it I ! | i M i M M ! I I I I | 1 i IN i ! I i i i i I i i M i M M !! i M i * M I * i i i M i i !! i 

m GENERATE OUTPUT DATA FILES n 

M i i l l ! i I 1 i i i < ! I < I ! i i l ! < 1 • i M i M i < l M ! ! ! I ! * i l ! ! ! M ! I i M < i I ■ l M i i * I < 

l 

PRINTER IS 7® 1 

i 



PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
FOR I 



USING RUN NUMBER" " , 1 5X , DO" ; Try 

USING POUER IN UATTS"" . 1 1 X , Z . D" ; Power 

USING SPACING IN MM" ", 1 2X , 1 ®A" ; Space* 

USING BLOCK #’ s POUERED" " . 8X , 1 ®A" ; B1 np$ 

USING AMBIENT TEMP. DEG. C" " , 5X , DD. DO" i T inf 

USING FILM TEMPERATURE"*' ,9X ,DDD. DO" ; Tf ilm 

USING FILM TEMPERATURE, C" " , 6X , DD. DD" i Tf 1 1 me 

USING THERMAL COND. OF WATER" '* , 3X , Z . DOD" i Kh2 o 

USING EXPANSION COEFFICIENT , B"** , 2X , 3DE" i Beta 

USING ""“KINEMATIC VI SCOSI TY“ " , 6X , 3DE“ ; K l nvi sc 
USING PRANDTL NUMBER" ", 1 1 X , DD. DO" •, Pr 



"BLOCK #" ."QCOND" ."QCONV" ,"GR 4" t “T*" , “H" , “NU 4" 
1 TO 8 



NU* 



1 18 



2010 Block-9-I 

2020 PRINT USING " 3X, D, 3X , S3DE , 3X , SZ . DOD , 3X , S3DE , 2X , S3D£ , 3X , S3DE , 2X , S3DE , 2X , S3 
DE“ » Block , Qcond( I ) , Qconv< I ) ,6ratar< I ) .Tatar < I),H(I),Nu(I) ,Nustar( I ) 



2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
21 10 
2120 
2130 
2140 
2150 
21G0 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
Dt f < I ) . 
2250 
2260 
2270 
2280 
2290 
2 300 
2310 
2320 
2330 
2340 
2350 
2360 
2 370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
Z500 



NEXT I 
PRINT M M 
PRINT “ ” 

PRINT " " 

PRINT •' " 

I ! ! m | I | i t I i M i i m ! ! ! M I ! M I m l I m M ! ! I M l l I I I l ! ! ! ! ! M !! ! M i M l m l l M 

!! CALCULATE TEMPERATURE DIFFERENCES M 

t l M ! | | M ! ! m I l | | M M m it I j | I u (I l ! ! j ! | i | l l l l i M ! ! M ! M ! ! ! i M M ; I { l m 

K *0 

FOR N-1 TO 8 
Dtf < N)*Tt< K ) 

Dtr<N)*Tt(K+1 ) 

Dtl<N)-Tt(K+2 ) 

Dtt(N)»Tt(K+3) 

Dtb< N )*Tt( K+4 ) 

Dth<N>«Tt(K+5) 

K-K + 6 
NEXT N 

PRINT “BLOCK 4“ , “OTF" , “DTR" , “DTL" , “DTT" , “DTB“ , "DTH“ 

FOR 1-1 TO 8 
Block *9-1 

PRINT USING “ 3X , D, 3X , S3DE . 3X , S3DE , 3X , S3DE , 3X , S3DE , 2X , S3DE , 2X , S3DE“ ; B1 ock , 
Dir< I ) , D*t 1 ( I ) ,Dtt( I ) ,Dtb( I ) ,Dth( I ) 

NEXT I 

M M ! i i I i | it I i m m I M t i | ! i I I M | ill I | i I M I I M i m m i i | j ! t M | i i I ii i I m I 

" CALCULATE NON DIMENSIONAL TEMPERATURE DIFFERENCES M 

i l i i l l M l i i l l !! i i l l ! i ll !! I i i l ! I M l ; l !| i I I l i ii i I l * l i M i * * M M ! M M M i 

FOR 1*1 TO 8 
! IF 1-1 THEN 
! Dthn< I )-88888 

! ELSE 

! IF 1-6 THEN 
! Dtf n< I ) *88888 

Difn< I )-Dtf < I )/Tstar< I ) 

Dtrn( I )-Dtr< I )/Tstar( I ) 

Dtln<I)*Dtl<I)/Tstar<I) 

Dttn( I )-D tt< I )/Tstar( I ) 

Dtbn( I ) -Dtb< I )/Tatar( I ) 

Dthn( I >-Dth< I ) /Tatar ( I ) 

IF 1-1 THEN 
Dt hn( I )-88888 
ELSE 

IF 1-6 THEN 
D*tf n( I )*88888 
END IF 
END IF 

NEXT I 
PRINT “ " 

PRINT “ “ 



2510 PRINT ”X/L” . "DTFN” . M DTRN” .“DTLN” . " DTTN” , “DTBN" , ''DTHN'' 

Z5Z0 Xs-.795 
Z530 6ap-.0ZS4/.3048 

Z540 FOR 1-1 TO 8 
Z550 Xl-Xs-< ( 1-1 )*G ap) 

Z5S0 PRINT USING " . DOO , 3X , S3DE , 3X , S3DE , 3X , S3DE , 3X , S3DE , ZX , S3DE , ZX , S3DE” ; X 1 , Dt f 
n( I ) . Dtrn( I ) ,Dtln( I ) # Dttn( I ) ,Dtbn< I ) ,Dthn< I ) 

Z570 NEXT I 
Z580 PRINT ” M 

Z590 PRINT " " 

ZG00 PRINT **DTR - RIGHT T. C. - TINF, etc.” 

Z610 PRINT ”T* IS TEMPERATURE SCALING FACTOR* 1 

Z6Z0 PRINT "GR U IS MODIFIED GRASHOF NUMBER** 

ZB30 PRINT **H IS HEAT TRANSFER COEFFICIENT” 

2640 PRINT *’ NU U IS NUSSELT NUMBER” 

Z650 ! OUTPUT 6RoadiQcond< • ) ,Qconv< * ) ,Gr< * ) ,Ttt< * ) ,B1< * ) ,Dtf ( • ) ,Dtr< * ) , 

Dtl< • > ,Dtt< • > .Dtb< • > ,Dth< * > 

Z660 1 

2670 I! NOTE THIS PROGRAM IS USED IN CONJUNCTION UITH THE 

Z680 !! THE DATA AQUISITION PROGRAM AND ACCOUNTS FOR THE 

Z690 ! ! INOPERATIVE THERMOCOUPLES ON BLOCKS 3i8 

Z700 ! 

Z710 END 
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1 ® 

z® 

3 ® 

4 ® 

5 ® 

6 ® 

7 ® 

8 ® 

9 ® 

I®® 

1 1 ® 

12 ® 

1 3 ® 

14 ® 

15 ® 

16 ® 

17 ® 

18 ® 

19 ® 

20 ® 

21 ® 

22 ® 

2 3 ® 

24 ® 

2 5 ® 

26 ® 

27 ® 

2 8 ® 

29 ® 

30 ® 

31 ® 

32 ® 

33 ® 

34 ® 

35 ® 

36 ® 

37 ® 

38 ® 

39 ® 

40 ® 

41 ® 

42 ® 

4 3 ® 

44 ® 

45 ® 

46 ® 

47 ® 

48 ® 

49 ® 

50 ® 

51 ® 

52 ® 



THERMOCOUPLE VARIATION WITH TIME PROGRAM 



M I ! II I M ! M I | M I M M l | I M | | ll ! ! M I M I M I ! I !! | ! I ! ! M ! ! ! ! ! I I ! M ) ! 

M TEMPERATURE TRANSIENT RESPONSE PR06RAM M 

i j j i i | i | | | m i | i ii i | i | i i ll | | i | | | | i i | |i I i | | | i | i i | i i i i i i i | i || ! | i it | 

REAL TZ ( 575 ) , T3( 575 ) , T4< 575 ) , T5< 575 ) . T6< 575 ) , T7< 575 ) , T8( 575 ) 

REAL V2( 575 ) ,V3<575) , V4< 575 ) , V5< 575 ) , V6< 575 ) , V7< 575 ) ,V8<575) 

REAL G( 575 ) 

N-575 

PRINTER IS 701 

CREATE BOAT " TRANS : , 7®0 , ® , 1 M , 5®Z 

ASSIGN STraill TO " TRANS : .70®,®, 1" 

| l i I I I i i i | I i ! i I ll i I i I I I i i i I ! I M I i i i i * I ‘I M t i 1 I M i I * 1 I i I I i I I M M I ! I I I I I I M I i 

BEEP 

FOR J-l TO 5 

OUTPUT 7®9 ; "CONFMEAS DC V, .USE 0" 

ENTER 7®9 { V 1 ( J ) 

OUTPUT 7®9 i "CONFME AS DC V .USE ®" 

ENTER 7®9 ; V2 ( J ) 

OUTPUT 709 {"CONFMEAS DCV. .USE 0" 

ENTER 7®9 ; V3( J ) 

OUTPUT 7®9 ; "CONFMEAS DCV .USE 0" 

ENTER 7®9 ; V4( J ) 

OUTPUT 7®9 ? "CONFMEAS DCV. .USE ®" 

ENTER 7®9;V5(J) 

OUTPUT 7®9 ; "CONFMEAS DCV. .USE 0" 

ENTER 7®9;V6(J) 

OUTPUT 7®9 {"CONFMEAS DCV. .USE ®" 

ENTER 7®9;V7(J) — 

OUTPUT 7®9 {"CONFMEAS DCV .USE ®" 

ENTER 7®9?V8(J) 

NEXT J 
BEEP 

FOR J-6 TO N 

OUTPUT 7®9 {"CONFMEAS DCV. USE ®" 

ENTER 7®9 ; V 1 ( J ) 

OUTPUT 7®9; "CONFMEAS DCV, USE ®" 

ENTER 7®9 ; V 2 < J ) 

OUTPUT 7®9; "CONFMEAS DCV. USE ®" 

ENTER 7®9 ; V 3( J ) 

OUTPUT 7®9; "CONFMEAS DCV, .USE ®" 

ENTER 7®9 ; V4( J ) 

OUTPUT 7®9{ "CONFMEAS DCV USE ®" 

ENTER 7®9 ; V5( J ) 

OUTPUT 7®9{ "CONFMEAS DCV USE ®" 

ENTER 7®9 ; V6( J ) 

OUTPUT 7®9 {"CONFMEAS DCV .USE ®" 

ENTER 7®9 ; V7( J ) 

OUTPUT 7®9; "CONFMEAS DCV .USE ®" 

ENTER 7®9;V8( J) 

IF J*N THEN 
BEEP 245®. 97* Tone ,1.5 
END IF 
NEXT J 
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530 A*. 0006797 

540 B=ZS8ZS. 1328 

550 0-607789.2467 

560 D--Z19SZ034. 3364 

570 E = 83708 10996. 1874 

580 FOR J«1 TO N 

590 IF J<5 THEN 

600 G( J ) *0 

610 ELSE 

6Z0 G( J ) *( J-5 ) *T in 

630 i Tim IS TIME BASED ON TIME PROGRAM RUNS FROM 2ND BEEP TO TONE DIVOED 

640 ! BY N-5, E.G. FOR 7 THERMOCOUPLES Tin IS .745, FOR 8 IT IS .852 

650 END IF 
660 NEXT J 

670 PRINT " " 

680 PRINT “ M 

690 FOR J-1 TO N 

700 T1( J)-A+(B»VI( J) ) + <C*< VI ( jrZ) )+(D*( VI ( J) A 3) ) + (£*( VI ( Jr 4) ) 

710 TZ< J)-A+< B*VZ( J ) ) + ( C*< V2( Jr 2 ) ) + ( D*< VZ< Jr 3) > + < E*< VZ < J) A 4) ) 

720 T3( J ) = A+( B*V3< J) )+<C*< V3< jrZ> )+( DM V3< Jr 3) )+(E*( V3( Jr 4) ) 

730 T4< J)=A+( B»V4( J ) ) + (C*( V4< J) A Z ) )♦< 0* < V4( J r 3) ) + < E*< V4< J) "4) ) 

740 T5( J ) = A+( B» V5< J ) )♦( C* < V5( J r 2 > )+(D»(V5(J)"3) ) + (E*(V5(J)"4) ) 

750 T6< J)«A+< B*V6< J) )+< C*< VG< J) A Z ) ) + ( D*< V6< J )*3) )+( E*< V6< J) A 4) ) 

760 T7< J)-A+( B*V7( J) )♦< C*< V7( J) A 2 ) ) + ( D*( V7( Jr 3> )+( E*( V7< J) A 4 ) ) 

770 T8< J)-A+( B»V8( J) ) + ( 0< V8< J) A 2 ) ) + < D*< V8< J) A 3) ) + < E*< V8< J) A 4) ) 

780 PRINT USING “IX, DDD , 1 X ,4DE , 1 X ,4DE , 1 X ,4DE , 1 X , 4DE ,1X, 4DE , 1 X ,4DE , 1 X ,4DE , 1 X , 4 
DE" ;G< J) ,T1( J) ,T2( J) ,T3< J) ,T4( J) ,T5(J),T6( J) ,T7( J) ,T8( J) 

790 NEXT J 

800 OUTPUT ©Trail 1 ;G< O ,T! < * > ,TZ< *> f T3< • > t T4< * > ,T5< • > ,T6< • > ,T7< • > .T8< • > 

810 END 
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APPENDIX B 



EQUATIONS FOR DETERMINING FLUID PROPERTIES 

p = -7.945E-10 T 2 FaM + 5.736E-5 T FaM - 9.781E-3 

v = 4.697E-9 T FaM - 2.537 E-6 T FllM + 1.342E-3 
PL = 3.235 E-7 T 2 FaM - 2147 E-3 T FILM + 0.03617 
v = VfJ. 

k = 1. 42 E-3T FaM +0.187 
Pr = 4.66E-3 T FaM - 2.92 T faM + 463 

Curve fits were determined for range of temperatures from 280 to 
300 K. 
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APPENDIX C 



SAMPLE CALCULATIONS 



The calculations were based on temperatures derived from a run 
with all blocks powered to 1.0 watt, no shroud in place. When 
appropriate, values are calculated for block number 6. 



A. CHARACTERISTIC DIMENSIONS 



Perimeter = 



'2(0.00787 + 0.02388) + (2(2(0.00787 + 0.00610)))' 
.+ (2(2(0.02388 + 0.00610))) 



= 0.2393 m 



/l, = (0.02388X0.00787) = L879 x 10 m 2 
A 2 = (0.02388X0.00610) = L457 x 10' 4 m 2 
A 3 = (0.00787X 0.00610) = 4.80 x 10" 5 m 2 
A™,* (1.879 x 10" 4 ) + (2)(l.457x 10" 4 ) + (2)(4.80 x 10' 5 ) = 5.75 x 10"V 

/ (J IAL 

A total* = (5. 75 X 10' 4 - L 879 x 10' 4 ) = 3.87 x 10" 4 m 2 

f~ 5. 75 X 10 1 rT ^ t 

L - 0. 2393“- 2 40 x 10 m 
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B. SURFACE AVERAGED TEMPERATURES 



(1 9.82)( 1.879 x 10" 4 ) +(19.48)(4.80x 10' 5 ) 

+ (19.49)( 4. 80 x 10 -5 )+ (20.62) (L457 x 10" 4 ) 
+ (20.26)( L457 x 10" 4 ) 



■ AVG ( 8 ) 



(5.75 x 10" 4 ) 



- = 20.08^ 



( 19.4 1)( 1.879 x 10" 4 ) + (18.99)(4. 80 x 10" 5 ) 

+ (19.05)( 4. 80 x 10" 5 )+ (19.88)(L457 x 10" 4 ) 
+ (20.08)(L457 x 10" 4 ) 



AVG ( 7 ) 



(5.75 x 10 -4 ) 



- = 19.64 °C 



(19. 16)( 1.879 x 10" 4 ) + (1 8.67)( 4. 80 x 10" 5 ) 

+ (18.6 1)(4. 80 x 10" 5 )+ (20.00)( L457 x 10" 4 ) 
+ (19.57)( L457 x 10" 4 ) 



AVG ( 6 ) 



(5.75 x 10" 4 ) 



= 19.40 °C 



(1 8.96)( 1.879 x 10" 4 ) + (18.33)U 80 x 10 -5 ) 

+ ( 1 8.47)( 4. 80 x 10" 5 )+ (19.68)(l457 x 10' 4 ) 
+ (19.37)(L457 x 10" 4 ) 



AVG ( 5 ) 



(5.75 x 10~ 4 ) 



- = 19.16 °C 



( 18.65)( 1.879 x 10" 4 ) + (18.30)(4. 80 x 10" 5 ) 

+ (1 8.39)( 4. 80 x 10" 5 ) + (19.75)( L457 x 10" 4 ) 
+ (I9.06)(l457 x 10" 4 ) 



AVG ( 4 ) 



(5.75 x 10" 4 ) 



= 18.99 °C 
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AVG (3) 



’(17.76)(4.80x 10' 5 ) + (1 8. 10)( 4. 80 x 10‘ 5 ) 

_+ (19. 1 8)( 1 457 x 10~ 4 ) + (18.87)( L457 x 10~ 4 )_ 
( 3.87 x 10' 4 ) 



= 18.77 °C 



( 18.07)( 1.879 x 10“ 4 ) + (17.40)( 4. 80 x 10" 5 ) 

+ (17.79)( 4. 80 x 10‘ 5 ) + (18.85)(L457 x 10" 4 ) 
+ (1 8.56)( L 457 x 10" 4 ) 



AVG (2) 



(5.75 x 10' 4 ) 



= 18.32 °C 



AVG (1) 



(17.73)( 1.879 x 10" 4 ) +(17.25)(4.80x 10" 5 ) 

+ (17.25)( 4. 80 x 10" 5 )+ (18.60)(L457 x 10" 4 ) 

,+ (1 7.78)( L 457 x 10~ 4 ) 

(5.75 x 10" 4 ) 



- = 17.89 °C 



C. TFILM 



f _ 20.08+ 19.64+ 19.40+ 19.16+ 18.99+ 18.77+ 18.32+ 17.89 = J9 



T 19.03+ 13.23 ir rT 

1 FILM 2 — IO. iJ ^ 



D. TEMPERATURE EXCESSES 

Using data from block 6, temperature excess is evaluated as: 

A7> = 19.16- 13.23 = 5.93°C 



AT r = 18.67 - 13.23 = 5.44°C 
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A T l = 18.61- 13.23= 5.38^ 



AT t = 20.00- 13.23 = 6.77°C 



A T„ = 19.57 - 13. 23 = 6.34° C 

D 



A T h =24.20- 13.23= 10.97°C 



E. CONVECTED HEAT FLUX 



_// 0.984 

q = — = 5240 — - 

L 879x 10 m 2 



F. TEMPERATURE SCALING FACTOR 



(5240)(z40 x 10~ 3 ) 

1 ~ (0.597) " ZL1 ^ 



G. MODIFIED GRASHOF NUMBER 



Gr*= 



(9. 8 1) ( 1 63 x 10~ 6 )(5240)(l40 x 10~ 3 ) 
(0. 597)(l 12 x 10' 8 ) 



H. HEAT TRANSFER COEFFICIENT 



*= = 849-* 



(19.40- 13.23) m 2 K 



= 371 
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I. NUSSELT NUMBER 



Nu 



(849)(2.40x 10~ 3 ) 
(0.597) 



= 3.41 
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